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Abstract. We propose, analyze, and demonstrate a discontinuous Galerkin 
method for fractal conservation laws. Various stability estimates are estab- 
lished along with error estimates for regular solutions of linear equations. 
Moreover, in the nonlinear case and whenever piecewise constant elements 
are utilized, we prove a rate of convergence toward the unique entropy solu- 
tion. We present numerical results for different types of solutions of linear and 
nonlinear fractal conservation laws. 



1. Introduction 
We consider the fractional (also called fractal) conservation law 

dtu{x,t) + d^fiu{x,t)) ^ gx[u{x,t)] {x,t) e :=Mx (0,T), 



^^'^^ "I u{x,0) = uo{x) xd 

where / is a Lipschitz continuous function and gx is the nonlocal fractional Laplace 
operator —(—9^)'^/^ for some A € (0, 1). This operator can be formally defined by 
Fourier transform as 

(1-2) 5Al^)](0 = -|eiV(0 

or, equivalently, by a singular integral (cf. [211 HZ]) as 

r / M f (p{x + z)- fix) 

gxmx)\ = cx / piYT^^ dz 

J\z\>0 V-\ 

for some cx > 0. For sake of brevity, we often write g instead of gx in the following. 

Nonlocal partial differential equations appear in different areas of engineering 
and sciences. For example, the linear nonlocal partial differential equation 

(1.3) dtu — d^u — dxU + u = gx[u] 

is a nonlocal generalizations of the famous Black-Scholes' equation in finance jl6j . 
and has received a lot of attention in the last decade. In recent years, attention has 
also been given to nonlinear nonlocal equations like 

(1.4) dtu + udxU ^ gx[u\, 



known as the fractional Burgers' equation. Equation (1.4) finds application in cer- 
tain models of detonation of gases (cf. |5D]) characterized by an anomalous diffusive 
behavior which can be described by means of the fractional Laplacian. We refer 
the reader to [21 El HI] , and the references therein, for further applications in hy- 
drodynamics, molecular biology, semiconductor growth and dislocation dynamics. 
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Many authors, see [2 13 El El H 13 [HI [21], have contributed to settle issues like 
well-posedness and regularity of solutions for the fractional conservation law 
In the case A € (1, 2), ( |1.1| ) is the natural nonlocal generalization of the viscous con- 
servation law dfU + dxf{u) — d^u. Such equations turn a merely bounded initial 
datum into a unique stable smooth solution (cf. [50]). The case A € (0, 1) is more 
delicate. Alibaud's entropy formulation is needed to guarantee well-posedness [2], 
and the solutions may develop shocks in finite time [3j; the diffusion is no longer 
strong enough to counterbalance the convection, and equation (1.1 1 fails to regu- 
larize the initial datum. In the critical case A = 1, Alibaud's entropy formulation 
is still needed to ensure well-posedness, however, solutions should be smooth as in 
the case A e (1,2) - see Kiselev et al. [25] for the case of the fractional Burgers' 
equation. 

A vast literature is available on numerical methods for nonlocal linear equations 
like ([O]). The interested reader could see, for example, [i] [TU l [TTl [T^l [T ^ [TS] [25] . 
However, numerical methods for nonlocal nonlinear equations like (1.1 ) are far from 
being abundant. Dedner et al. introduced in [17j a general class of differences meth- 
ods for a nonlinear nonlocal equation similar to ( 1.1 1 coming from a specific problem 
in radiative hydrodynamics. Droniou |19j was the first to analyze a general class 
of difference methods for (1.1), he proved convergence toward Alibaud's entropy 
solution, but produced no results regarding the rate of convergence of his methods. 

In this paper we study a discontinuous Galerkin (DG) approximation of (1.1). 
The DG method is a well established numerical method for the pure conservation 
law dtu + dxf{u) — 0. Some of the important features of this method are stability 
and high-order accuracy. Moreover, when piecewise constant elements are used, 
the DG method reduces to a conservative monotone difference method (cf. [25] ) 
which converges to the entropy solution with rate 1/2 (cf. the well known results 
of Kuznetsov 02] )■ For a detailed presentation of the DG method for pure conser- 
vation laws, we refer to Cockburn [14 . 

In this paper we propose a DG approximation of (1.1) in the case A G (0,1), 
and prove that we retain the main features of the DG method in our nonlocal 
setting. We show L^-stability, and prove high-order accuracy for linear equations. 
Moreover, when piecewise constant elements are used, we derive two fully discrete 
numerical methods, an implicit-explicit method as in and a fully explicit one, 
and prove convergence toward a BV entropy solution of ( |1.1[ ) (cf. Definition 4.1 



below) with a certain rate. For the implicit-explicit method we prove convergence 
with rate 1 /2 while for the fully explicit one we prove convergence with a lower rate, 
min{l/2, 1 — A}. To prove the rate of convergence, we generalize the Kuznetsov 
argument |26j to our nonlocal setting, and, as a byproduct, we obtain the following 
theoretical result: Alibaud's entropy formulation and the BV entropy formulation 
are equivalent whenever the initial datum is integrable and of bounded variation. 

Finally, several numerical experiments have been performed to illustrate the de- 
veloped theory. Among other things, we are able to reproduce the theoretical results 
(absence of smoothing effect due to persistence of discontinuities and formations of 
shocks) obtained in [3 [25] for the fractional Burgers' equation. 



2. A SEMIDISCRETE DG METHOD 

Let us introduce the space grid xi = iAx, i e Z, and let us label li = (x^jXi+i). 
We call P^{Ii) the set of polynomials of degree at most k e {0, 1,2,.. .} with support 
on the interval and consider the Legendre polynomials (cf. [14] for details) 

WQ,i,'Pis, - ■ ■ iVk,i\i G P'\Ii) for aU j = 0,...,fc. 
Each ip S P^[Ii) is a linear combination of the functions {</3o,i, • • ■ , 'Pk,i}- 
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If we multiply (1.1) by an arbitrary </? £ P^{Ii)^ integrate over the interval 
integrate by parts, and replace the flux / by a numerical flux F, we get 



(2.1) 



fiu)ip^, + Fiu^+i)(p{x^^-^) - F{ui)ip{xj) = / g[u]ip 



As usual for DG methods, the numerical flux F{ui) — F{u{x^ ), u{x^)) satisfies the 
following assumptions: 



Al 
A2 
A3 
A4 



F is Lipschitz continuous on M x M, 
F{a, a) = f{a) for all a G M, 

F is non-decreasing with respect to its first variable, 
F is non-increasing with respect to its second variable. 



The goal is to find a function u : M x [0, T] ^ M, 



(2.2) 



H2;,i) = ^^UpA{t)(ppAx), 



which satisfies (2.1) for all ip G P''{Ii), i G Z. Let us fix ip{x) = '^„^o<^q,i'Pq.i{x), 



and plug (2.2) into (2.1) to get 

k 



?=0 



Ax d 



where F{ui) — -F'(X]p=o Sp=o derive the above expression 

we have used some well known properties of the Legendre polynomials: for all i G Z, 

ioT p — q 



^p^i^q^i 



dx 



otherwise ^ ^ ^ 



where we have denoted with (/?(a:^), (p(x^) the (right and left) limits of (/^(s) as 
s ^ Xi. The semidiscrete method (i.e., discrete in space and continuous in time) 
we study is the following: for all g = 0, . . . , fc and i G Z, 



(2.3) 



UqA^) = ^ Jj.Uo{x)^q4x)dx. 



3. Nonlinear L^-stability and convergence in the linear case 

Let V'' := {u : G P''{Ii) for all i G Z} be the space of piecewise polynomials, 
and let 7J^/^(R) be the fractional Sobolev space with norm 



2 



II l|2 — II l|2 I |2 ^ I |2 — f f W{z) — u{x)] 

ll"ll_ff^/2(R) II"IIl2(r) + pl_H-A/2(R) and |m|^a/2(h) •— / / u _ ^n+A dzdx. 

JRJR 1^ -^1 

Let us note that the space H'^^^(R) also contains discontinuous functions (cf. [21 
Lemma 6.5]). Moreover, let us denote with H^^/^{R) the dual space of iJ'^/^(R 



and let us point out that, as shown in the proof of Corollary A. 3 below, g[u] G 
iJ~'^/^(IR) whenever u G ff'^/^(M). In the following, all the integrals of the form 
J-^g[u]v, where the functions u,v £ i7'^/^(M), should be interpreted as the pairing 
{g[u\,v) between H^^^{R) and its dual. 
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Theorem 3.1. (Stability) If also /(O) — 0, then any solution u of (2.3 I belonging 
to Ci([0,r];i7-^/2(M)) IS L"^ -stable: 

l|w(->7')lli2(R) +CA / |M(-,OI?fA/2(R)rfi < ||uo|li2(R). 







The above result generalizes a well known result for the DG method for pure 
conservation laws (cf. |T3| Proposition 2.1 and Theorem 4.2] for details). 



Proof. By construction, u{-,t) satisfies (2.1) for all test functions if G P {li). Let 



us choose the test function <p = u{-,t), sum over i € Z, rearrange the terms in the 
sum and integrate over time to get 

/ utii^ V'LF(Mi)({i(a;+)-{t(a;,^))+ / f{u)ux + / 9[u]u. 
Jqt Jo •'^i •'Qt 

Due to the assumptions made (including A1-A4), each term in the above expression 
is well defined. The first term is clear while the last term is well defined by Corollary 



A. 3 The remaining terms makes sense for all functions in V'' L^(M) since the 
point values are well defined. To see this, note that, since /(O) = and / is 
Lipschitz continuous, f{u) and F{u) belongs to L^(IR) since u does. We can then 
conclude, using the Cauchy-Schwarz inequality, if the function v, v = in Uig^-^i: 
belongs to L^(M) (this is the regular part of the distribution u^)- But this again is 
an easy consequence of the regularity of the Legendre polynomials ipp^i and their 
othogonality which implies that 

^^Cp,i(0<^p,j(a;) e L^{Qt) if and only if / 4A*)^* < °°- 

iez p=o iezp=0''0 

Let us now prove stability. Since 

fu(x) /'u(x^^^) r'^i^t) 



g[u\u. 

Qt 



^ f{u)Ux = J^iJ f)x= j / - / /, 

we find that 

[ UtiL ^ \nn^mx+) - U{xi)) - r^^^ ^ f{x)dx] + f 

JQt jo Ji~<^,) JQ 

It is well known that a fiux satisfying A2-A4 is an E-flux (cf. [H]), i.e. 

F{u,){u{x^) - u{x:r)) - / f{x)dx < for aU i e Z. 

Thus, by Corollary [X3l 

2M-,T)\\l2(^m} + Y J^J^ |^_^|i+A dzdxdt< -WuoWmR), 

and the proof is complete. □ 

In the linear case, equation ( |1.1[ ) reduces to 
(3.1) dtu + cdxU = g[u] 

where c € M. Let us prove the following result. 

Proposition 3.2. Let uq e H''^^{M.). k > Q. Then, there exists a unique function 
u G H^'^^{Qt) which solves (3.1). Moreover, 



(3.2) ||u(-,<)||^fc+i(R) < \\uo\\ifk+i( 
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Proof. Since (3.1 ) is linear, its Fourier transform, dtu + i^cii — ~\^\^u, has solution 

This implies existence plus, using Plancherel theorem, L^-stability and uniqueness. 
L^-stability for (weak) higher derivatives can be obtained as follows: take the de- 
rivative of (3.1), repeat the above procedure, and iterate until the k-th derivative. 
Regularity in time can be shown by using equation (3.1) and regularity in space. □ 

As pointed out by Cockburn |14j . in the linear case all relevant numerical fluxes 
(Godunov, Engquist-Osher, Lax-Friedrichs, etc.) reduce to 



^{a + b)-f{b^a). 



(3.3) F{a,b) 

We use this flux to prove the following result: the order of the semidiscrete method 
(2.3) increases along with the degree k of the polynomial basis used. 



Theorem 3.3. (Convergence) Let u £ H^^^{Qt), k > 0, be a solution of (3.1) 
and u G C^([0, T]; i7'^/^(M)) be a solution of the semidiscrete method (2.3). Then, 
there exists a constant Ck,T > such that 

||M(-,r)-?i(-,r)|U2(R) <cfc,TAx'=+i 

The above result, called high-order accuracy, generalizes a well known feature of 
the DG method for pure conservation laws (cf. |14i Theorem 2.1]). We are able to 
prove this result since, as shown in the proof below, the error due to the local terms 
(cfe^TAa;*''+^/^) is bigger than the one due to the nonlocal term {ck^T^x'''^^~'^^'^)- 

Proof. By construction, for all test functions ip E V'^ D -L^(M), 



CUipx 



g[u\ip. 



Note that u satisfies the analogous expression 
(3.4) 



/ utip + Y\F{ui){p(x^) - ip{xl)) - / cuip^ = / g[u]ip. 

JR L J li J Jr 

To prove the above relation, let us multiply (3.1 ) by a test function tp and integrate 
over li. Note that, thanks to the iJ'^+^-regularity of u, u is continuous (by Sobolev 
embedding). Thus, since F satisfies assumption A2, we get that 



Ut^p -f cuj,p - g[u]p 



Ufip- / cuip^ + F{u^+i)ip{x^.-^) ~ F{ui)(p{xf) 



g[u]ip. 



We obtain (3.4) by summing over all i £ Z and rearranging the terms in the sum. 



Let us introduce the bilinear form 

B{e,p) 



et(p + 

where e := u ~ u £ H^''^{R). Let us call u the L^-projection of u into V'': i.e 
(u(a;) — u{x))(pji{x)dx = for all j = 0, . . . ,k and i G Z. 
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Note that, by Lemma |A.4[ u € n L'^{R) implies u e H^^'^(R). Let us call 
e u - ?i e H^^'^iR). Since B{e, e) = 0, B{e, e) = B{e - e, e) = B{u - u, e) or 

[ [ ete^ [ [{u-u)te-[ V [F(e,)(e(a;7) - e(a;+)) - /" cee^] 

Jo Jo JR Jo '- J li 



[^((u - ").)(e(a:,r) - e(a;+)) - / c(u - u)e, 

+ / / .9[e]e- / / g[e~e]e. 
Jo Jr Jo Jm 

Note that, since both e,e e H^/^{R), each term in the above expression is well 
defined (cf. the discussion in the proof of Theorem 3.3). One can argue as in [HI 
Theorem 2.1] to bound the local terms by c^.^Ax^'"'"''"'". Hence, 

[ [ ete< Cfc,TAa;2'=+i + [ [ g[e]e - [ [ g[e - e]e. 
Jo Jr Jo Js. Jo Jr 

Let us denote by I wha t it i s left to estimate on the right-hand side of the above 
inequality. By Corollary A. 3 the i/'^/^-regularity of both e, e implies that 



1 = 



< 



T 

JM 



9[e]e+l: [ [ 9[e]e-l^ [ [ g[e ^ e]{e ~ e) 
Jo Jr ^ Jo Jr 



and, by Lemma | A. 5 [ 



(u- u)(-,t)||^V2(R) < Cfc||M(-,t)||^fc+i(R)A; 



2fe+2-A 



Thus, using the i/'^+i-stability of u, /J'/gete < Cfe,T[Aa;2'=+i + Ax^'^+^-^J, and, 
since e(a;,0) — and ||e|| = — u) — (u — u)|| > ||e|| — ||w — u||. 



l|e(-,r)|||2(K) < Cfc,T 



Ax^'^+i + Ax2'=+2-^ + Aa;2^-+2 < c^.tAo^^'^+i. 



□ 



Remark 3.4. Let us prove that a solution u G (^^([0, i]; i7'^/^(M)) of the semidiscrete 
method (2.3) actually exists up to some time t > 0. We consider the map 



u{-,t) e n i^(M) ^ J^l' (i) the right-hand side of (2.3|, 



and call J^u{-,t) :— X^iez Eg=o Note that, using CoroUary A. 6 (here 
the assumption /(O) =0 is needed), 

(3.5) m(-, t)ev''n L^{R) => Tu{-,t) ev'^n L^iR), 

and, since both (/, F) are Lipschitz continuous, there exists a constant c > such 
that, for all u,v€V'' € L^iR), 



(3.6) 



\\{Tu - J^v){-,t)\\L2(R) < CIK-U- w)(-,<)||i2(R). 



Therefore, thanks to ( |3.5[ ) and (3.6), an application of the Cauchy-Lipschitz's the- 
orem yields the existence of a time t > and a unique solution 

ueC^{[0,t];V''nL^{R)) 



of the semidiscrete method (IZS]). To conclude, note that V'' n L'^{R) C H^/'^(R) 



by Lemma A. 4 
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4. Convergence in the nonlinear case 

We study the nonlinear case by using only piecewise constant elements (fc = 0): 

{fo,i, ■ ■ ■ , fk.i} = Woa}^ VO,i = 1/., 
where 1/. : M — > M is the indicator function of the interval li = {xi,Xi+i). Start- 
ing from the semidiscrete method (2.3), we derive two fully discrete methods: an 
implicit-explicit method and a fully explicit one. By adapting Kuznetsov's tech- 
nique |26| to our nonlocal setting, we prove that both methods converge toward a 
BV entropy solution of with a certain rate (cf. Theorem 4.4). In Corollary 
|4.5[ we show how this result ensures well-posedness for BV entropy solutions of 
(1.1). Note that, in the nonlinear case, even when pure conservation laws are con- 
sidered, no results concerning the rate of convergence are available for high-order 
polynomials (fc > 0). 

Let us introduce the time grid i„ — nAt, where n = {0, . . . , N} and NAt — T. 
We discretize the semidiscrete method (2.3) in time to obtain the implicit-explicit 
method 

1 



(4.1) 



< 



U? 



and the fully explicit one 

jjn+l 



Ax 
1 

~ Ax 



AtD^F{u:^,UT+^) 
UQ{x)dx, 



Atg{U 



n+l\ 



AtD^F{ur,ur+^) 

UQ{x)dx. 



AtgiU"") 



Here we have introduce the shorthand notation D^F{Ul\ Ul\^) := -s^iFiUl", Ul\^^)- 
F{U"_i, J7")) and the nonlocal operator 

1 glUn^dx^±.J2G]U^, 



Ax 



where := g[li.]dx (we denote with [/' 
by the grid values {U^}i£i, such that U"-{x) 

Proposition 4.1. For all G Z x Z, 



M M the step function generated 
J7" for all X e [xi,Xi+i)). 



feGZ 



\Gl\ < oo. 



Gl - 0, 



Moreover, G* > whenever i ^ j , while 

G\ = ~d\Ax^^^, where d\ := cx | 

Proof. See the appendix. 

Let us introduce the CFL condition 
(4.3) (Fi + F2) 



(J J — Cr-, 




f dz 









> 0. 



□ 



At 
Ax 



< 1 



for the implicit-explicit method (4.1) (here Fi,F2 are the Lipschitz constants of F 
with respect to its first and second variable) and the CFL condition 



(4.4) 



(Fi + ^^2) 



At 



dx 



At 
Ax^ 



< 1 



for the fully explicit method (4.2). In what follows, the relevant CFL condition is 
always assumed to hold. 
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Let US introduce the time discretization into (2.2) as follows: 

(4.5) u{x,t) = C/" for all {x,t) e [xi,Xi+i) x [t„,t„+i). 

Theorem 4.2. Let uq e L^(IR) n BV{M.). Then, both the implicit-explicit method 
(4.1 ) and the fully explicit method (4.2 ) enjoy the following properties: for all t > 0, 

i) ll"(-,OIU~(K) < ||mo||l-(r), 

\\u{-,t)\\L^K) < ||wo||li(R), 
Hi) \u{-,t)\BV{R) < \uo\bv{r)- 
Moreover, there exists a constant c > (whose value is independent of the dis- 
cretization parameter Ax) such that, for all s,t >0, 

iv) \\u{-,s) - t2(-,<)||ii(R) < c(|s - t| + Ax). 



Proof. We give here the proof for the fully explicit method (4.2). The proof for the 
implicit-explicit method (4.1) can be found in the appendix. 

Let us point out two consequences of Proposition |4.1[ In the first place, note 
that the fully explicit method (4.2 1 is conservative. Indeed, since J^jez l^jl °^ 
for all i E Z, 



(4.6) 



whenever 



EEiG}^i-Eif^iEiG}i<oo, 

< oo. Thus, since J2iez ^) ~ ^ ^^^^ ^ ^' 



which implies X^ipz ^t"^^ — X^igz^j"- second place, note that the fully 



explicit method ( |4.2[ ) is monotone in view of the CFL condition (4.4). 

We are now ready to prove the theorem. Indeed, monotonicity and Proposition 



4.1 



G|, = 0) imply item i. The proofs of items ii and Hi follow, word by 
word, the ones in j23t Theorem 3.6]. Finally, note that, since the numerical flux F 
is Lipschitz continuous in both variables, there exists a constant c > such that 



n+l 



At 



= Ati?_F(;7r, c/;;i) + ^ E ^i^^ 



(4.7) 



iez 



Ax' * 



At I 
Ax I 



jez 



Let us multiply both sides of (4.7) by Ax, and sum over all i ^7^. Since 

ElE^i^ji^ / \9m\dx^cxc\m\i'ii^^ - 

< caC||uo|1^7(r)|woIw(: 



\BV( 



(cf. Lemma (A.l)), we get ||[/"+ — [/"||li(r) < cAt which implies iv via (4.4). □ 



Let us introduce the definition of BV entropy solutions of (1.1). Let rjkiu) := 
\u-k\, T]'^{u) sgn(u - k) and qk{u) r^'^{u){f{u) - f{k)). 



Definition 4.1. A function u G L°°{Qt) is a BV entropy solution of ( 1.1 ) provided 
that the following two conditions hold: 

i) ue C{[0,T];L\R))nL°°{0,T;BV{R)); 
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a) for all k CzR and all nonnegative ip £ C^{Qt), 
rikiu)^t + qk{u)Lpx + r]'^.{u)g[u](pdxdt 



(4.8) 



Qi 



rik(ua{x))ip{x,0)dx ~ / rik(u{x,T))ip{x,T)dx >Q. 
: JR 

The nonlocal term in the above definition is well defined since, by the regularity 



of u, g[u\ is integrable over the domain Qt (this is a consequence of Lemma A.l I. 
Note that sufficiently regular solutions of ( |1.1[ ) are solutions according to the above 
definition while solutions according to the above definition are weak solutions of 



( 1.1 ) (this can be easily proved by choosing k as the supremum of We refer the 



reader to Alibaud's paper P] for the precise definition of a weak solution of (1.1). 

As already mentioned in the introduction, Alibaud's entropy formulation en- 
sures well-posedness for all bounded initial data. We prove that the BV entropy 
formulation is well-posed for all initial data belonging to a smaller set, the set of all 
integrable functions of bounded variation, and, therefore, Alibaud's entropy formu- 
lation and the BV entropy formulation are equivalent whenever the initial datum 
lies in this smaller set. 

The following lemma generalizes to our nonlocal setting a result due to Kuznetsov 



[26], and it is used in the proof of Theorem 4.4 Let us introduce the function 
if{x,y,t,s) = uj^{x — y)ujs{t — s) where uj^ € C^(M), a > 0, can be built as follows: 
choose w G C^(M) such that < w < 1, uj{x) — for all > 1 and Jg^uj{x)dx = 1; 
finally, call ujaix) := uj{x/a)/a. 



Lemma 4.3. Let uq £ L^{R) n BV{R), u be a BV entropy solution of (1.1) and 
u : Qt — > M 6e any function such that items ii-iv in Theorem\4-.^ hold. Let 



A[u, (fi, k] := / r]k{u)(pt + qk{u)'Px + r]'^{u)g[u]ipdxdt 

rik{uo{x))ip{x,0)dx - / r]k{u{x,T))ip{x,T)dx 



Qi 



and A^^s[u,u] := Jq^ A[u,ip{-,y,-,s),u{y,s)]dyds. Then, there exists c > such 
that, for all e > and < d < T , 

\\u{-,T) - i2(-,T)||ii(R) <c{e + S + Ax) - A,j[u, u]. 

Proof. See the appendix. □ 

The above Kuznetsov type of lemma allow us to prove the following rates of 
convergence. 



Theorem 4.4. Let uq e L^{R) n BV(R) and u be a BV entropy solution of ( 1.1 1 



a) If a is the solution of the implicit- explicit method (4.1), then there exists a 
constant ct > .such that 

\\u{-,T) - u{-,T)\\Lum < ctVAx. 



b) If u is the solution of the fully explicit method (4.2 1, then there exists a 
constant ct > such that 



h(-,r)-w(.,r)|Ui(R) <ct(Va^ + a 



X 



The rate of convergence obtained for the implicit-explicit method (4.1 1 general- 
izes to our nonlocal setting the rate of convergence obtained by Kuznetsov in |26) 
for local difference methods for pure conservation laws. We suspect the convergence 



rate for the fully explicit method (4.11 to be suboptimal. Anyway, to the best of 



our knowledge, no convergence proof for the fully explicit case was available in the 
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literature up to now (cf. Droniou |19) for an alternative convergence proof, without 
convergence rate, for the implicit-explicit case). 



Proof. The plan is to estimate —A^^s[u,u], and, then, use Lemma 4.3 to conclude 



Proof for the implicit- explicit method. Let us introduce the notation a A b = 
min{a,6}, aV6 = max{a,&}, rjf = |C/f -u| and = /(J/f V w) - /(t/f A m), where 
u — u{y,s). Note that —A^^s[iL,u\ can be rewritten as 



-A, 



■'Qt [ „=o lez -^^^ 



(4.9) 



ip{x,tn+l)dx 



+ {<i7-€-i)j^ v{x..,t)dt 



ri'^{u)g[u]ipdxdt \dyds. 



Indeed, using summation by parts, 



E / / V?Vt{x,t) + q^ip^{x,t)dxdt 
+ rll ^{x,0)dx~ij^ 



^^E"! E / [ip{x,tn+i,) - ip{x,tn,)]dx 

[ip{Xi+i,t) - (p{Xi,t)]dt 

+ ^2 ^{x,0)dx-'qf ip{x,T)dx 

J Xi J Xi \ 

^EE hr+'-'^D / v{x.t^+i)dx 



ieZ (. n=0 •'■^i 

+ E^r/ 



+ {q1-qU) V{x,,t)dt 



i£l n=0 



Let us exploit monotonicity to get 

C/f+i V fc < C/f V fc - A<i?_^^(C/r V fc, V fc) + Aa(fe,+^)(C/f +i).g(C/"+i)., 
A fc > C/f A fc - MD^F{UT A fc, t/f+i A fc) + Aa(_^,,.) (f/f +i)5(C/"+i),. 

Let us call := F(;7f V fc, C/;!^i V fc) - i^(C/f A fc, J/^!^! A fc), and note that, since 
|a - 6| = a V 5 - a A 6, we can subtract J7"^^ A fc from U]^^^ V fc to obtain the cell 
entropy inequality 

(4.10) ^^"+1 - rj- + ^{Q- - QU) - Mi^',{U:+')g{U-+'), < 0. 
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If we plug the above inequality into (4.9), we find that 

-'Qt [ n= 

At 



r ( 



ip{xi,t)dt 



Aa 



N-l 



i 



n=0 i£Z 

ri'^[u)g[u]ipdxdt \dyds. 



V{x,tn+i)dx 



Next, the right-hand side of the above inequality needs to be estimated. To this 
end, let us point out that, as proved in |23', Example 3.14], 



X E E 



At 



<l"-i) I (p{x,,t)dt 

dyds < ct 



Ax Ax 

— + ^ 
e 



Let us call tp" := ^ J^"''^ vi^i tn)dx and ip the step function built from {</?"} by 
taking (p{x,t) = ip"^ for all G x [tn,tn+i). Moreover, let us call J the 

term which still needs to be estimated, 



J 



ri'^(u)g[u]Lpdxdt ^dyd. 



Since g{U^^)i = ^ /^.'^^ g[U^]dx, we can rewrite J as 

nT+At 



J = 



'Qt [J At 

which can be split into Ji + J2 ~ Js, where 



ri'^{u)g[u]ipdxdt ~ / r]'^{u)g[u]ipdxdt >dyd. 



I. 



Ji := 



Vuiu)g[u]{'P ~ ip)dxdt \dyds, 
Qt J 

T+At 

ri'^{u)g[u](pdxdt ^dyds, 



At 



r][^{u)g[u]Lpdxdt jdyds. 



By Lemma A.l and Theorem 4.2 g[u] e L^{Qt) and, thus, both 

nT+At 

J2 < 



[ \9m< 


/ (pdyds > 


Jr 


^ Jqt J 



At 



\9M 



ipdyds jdxdt 
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are of order Ax (here, as the the following, we use the CFL condition to pass from 
At to Aa;). Moreover, 



Ji < 



/ \9m< 


/ \V- v\dyds \ 


Jqt 


^ JQt J 



Ax Ax 



since, for all (x, t) e Qt: there exists a constant c > such that 

'Ax Ax' 

\Lp — Lp\dyds < c [ 



(4.11) 



We now prove ( |4.11 1. Let us call oo^ the step function built from {w^.i}, j 
Sx XT'^^ a;e(s)ds, as follows: We(x) = uj^^i for all x S First, note that 



(4.12) 
Indeed, 



\uje{x) - uj^{x)\dx < Ax\uj^\bv{R)- 



e{x) - UJe{x)\dx = ^ 

= E 



|a;c(x) — U!e{x)\dx 



1 

Ax 



dx 



< 



Ax 



E 



|We(s) — UJeix)\dsdx 



< AxIwelsyfB)- 



dsdx 



Next, we note that for all {x,t) £ Qx, 



(4.13) / \ip-^p\dyds^ / / \uj^{x - y)ujs{tn - s) ~ uJe{x - y)ujs{t - s)\dyds, 
Jqt Jo Jr 

where i„ is such that t e (i„,i„+i). Moreover, using (4.121, 
(4.14) 



l^eix " y) ~ uje{x - y)\dy = / \oj^{y) - uj^{y)\dy < Ax\uj^\ 
while, since t e (t„,i„+i). 



Ax 



\BV( 



(4.15) 



\^^s{tn - s) - ujs{t - s)\ds < At\ujs\BV( 



Thanks to the estimates (4.14) and (4.15), an application of the triangular inequal- 
ity to the right-hand side of (4.131 yields (4.11) 

The above estimates e 
use Lemma llTsl to obtain 



The above estimates ensure that —A^,s[u,u\ < ct(— + ^). Therefore, we can 



A.T Ax 

\u{-, T) - u{; T)IUi(R) < CT ( e + <5 + Ax + — + — 



The conclusion follows by setting e = (5 = V Ax. 

Proof for the fully explicit method. Let us exploit monotonicity to get 

C/f+i V fc < C/f V fc - AtD^F{U: V fc, V fc) -f Aa(fe^+oo)(C/f +i)g(C/")„ 

C/f +1 Ak>U^Ak- AtD^F{Ur A fc, Ul\^ A k) + Aa(_oo,fc) (^/^+')5(C^")^■ 
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Proceeding as done in the proof for the implicit-explicit method, we obtain the cell 
entropy inequality 

- V? + ^{Q'l QU) AW,{u-+')g{u-), < 0. 

Let us add and subtract Atr]'f.{U"^^)g{U'^~^^)i to the left-hand side of the above 
inequality, and let us use the fact that the operator g{-) is linear to obtain 



n+l n 



At 



AW,{U-+')g{U- - - At4{U-+')g{U"+'), < 



If we plug the above inequality into (|4.9|), we find that 



-A, 



At 
Ax 



iQ7 ~ QU) 



2:1+1 



N-l 



Lp{Xi^ t)dt 
ip{x,tn+l)dx 



n=Q i£Z 
N-l 



ip{x,tn+i)dx 



AtJ2Y.'^'JU-+')g{U--U-+'} 
r/^{u)g[u]Lpdxdt >dyds. 



Xi + l 



(p{x,tn+l)dx 



Qt 

The only term left to estimate is 



/ AtAx Y.'^'u{Ur^')9{U- - U-+%^{x,,tn+i) \dyds 

''Qt [ n=0 ieZ J 

< f lAtAxY,T.\9{U"-U^'-')Mx^,t„+Adyds 

J Qt I ^— n I 



< 



n=0 ieZ 

\g[u[x,t)-u{x,t + At)]\< 



tfdyds \dxdt. 



Note that, using Lemma A.l and Theorem 4.2 (item iv), the right-hand side of the 
above inequality is easily seen to be of order Ax^^^ . 



Finally, using Lemma [4~3 



|w(-,r) -u(-,T)||ii(K) <CT[e + 5 + Ax + Ax 



„1-A 



Ax Ax 



e S 

and the conclusion follows by setting e ^ 6 = ^/ Ax. □ 

We conclude this paper by proving the following result which is a consequence 
of Theorem 4.4 the definition of a BV entropy solution of (1.1 ) is well-posed. 

Corollary 4.5. Letu^ £ (R) n BV (R) . Then, there exists a unique B V entropy 
solution of 



Proof. Let us give the proof using the implicit-explicit method (4.1). Needless to 
say, the fully explicit method (4.2 1 would also do. 
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Uniqueness. Let us assume that both u and v are BV entropy solutions of (jOJ. 
If we add and subtract the solution of the implicit-explicit method (4.1 ), we obtain 

\\u{-,T) - w(-,T)||ii(R) < \\u{-,T) - ?i(-,r)||ii(R) + \\v{-,T) - {t(-,r)||ii(R) 



which, by Theorem 



4.4 



is less than or equal to ctv Ax + ctv Ax for all Ax > 0. 
Therefore, uniqueness follows. 

Existence. Using a standard argument (cf., for example, [231 Theorem 3.8]), 
Helly's theorem yields the existence of a subsequence m — >■ m in L\^^ {Qt) as Ax — >■ 0. 



Moreover, uG C([0, T]; Li(M))nL°°(0, T; BV^(E)) by Theorem 4.2 To prove that u 



satisfies the entropy inequality ( 4.8 1 , we start from the cell entropy inequality (4.10). 
Let us choose a nonnegative test function e C^{Qt) and call ipf := ip{xi,tn). If 
we multiply both sides of (4.10) by ip" > 0, sum over i and n, and use summations 
by parts, we find that 



n+l 



n=l i£l 



At 

N 



n=0 ieZ 



Ax 



Ax^{^Or;°-^f^f}>0. 



A standard argument shows that all the local terms in the above expression converge 
to the ones appearing in the inequality (4.8), cf. e.g. [231 Theorem 3.9]. Let us now 
consider the nonlocal term. Note that (here (p is as in the proof of Theorem 4.4) 

N 



n=0 ieZ 
N 



= AxAt ^ Vkiur^')9{u"+').{v? - ^r') + 

n=0 igZ 

where, since there exists a constant c > such that \ip" 



T+At 
At 



N 



< cAx 



ri[.{u)g[u\(pdxdt 

R 

"+^1 < cAxforall {i,n), 
\g[u{x,t + At)]\dxdt. 



AxAtJ2Y.'^k{Ul'^')9{U'^^'h{v2-^-i 

n=0 iGZ 

Since g[u] G L^{Qt), the right-hand side of the above expression is of order Ax. 
To conclude, we prove that there exists a subsequence {ii} such that 

rT+At 



(4.16) 



At 



Vk{u)9[u\'pdxdt 



Ax^O 



Vki'u-)9[u]'pdxdt 



for a.e. A: G M. This is a consequence of the dominated convergence theorem since 
the left hand side integrand converges pointwise a.e. to the right hand side inte- 
grand. Indeed, first note that (f ^ ip pointwise and that a subsequence u — >■ u 
a.e. in Qt- Moreover, for a.e. fc G M the measure of {{x,t) G Qt ■ u{x,t) — k} is 
null. This means that rj'/^iu) — )• i]'i^{u) a.e. in Qt, since t^J, is continuous on M\{A;}. 
Finally, by Theorem |4.4[ 

rT 



\g[ii — u\\dxdt < c 



l-A 



dt < cxAx 



for all Aa; > 0, and hence a subsequence g[u] — >■ g[u] a.e. in Qt. The proof for all 
fc G M follows the one given by Droniou in [T^], and this completes the proof. □ 
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■O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



(a) T = 0.5 



(b) T = 1.3 



-O.B -0,6 -0.4 -0,2 □ 0.2 0,4 0.6 O.B 



-0,B -0,6 -0.4 -0,2 □ 0.2 0,4 0.6 O.B 



(c) T = 0.5 



(d) T = 1.3 



Figure 1 . Initial data (piecewisc linear) and solutions of the pure 
fractional equation (A = 0.5) with fc = and Ax = 1/160. 



5. Numerical experiments 



We have implemented the numerical method (2.3) in the cases fc = 0, 1,2 with 



fully explicit time discretization. To perform computations, we have set our nu- 
merical solutions to zero outside the region = {[x,t) : \x\ < 3/2, t > 0}. In other 
words, we have computed the value Up^i{tn+i) using only the values {C^p,i(ira)}, 
where Xi G il and p — 0, . . . , k. This has been done also at the boundaries |x| = 3/2. 



Remark 5.1. Due to infinite speed of propagation (cf. [2 ), solutions of ( 1.1 1 do not 



have, in general, compact support. Therefore, the use of the region D, introduces 
an additional error which we have not considered in Theoreml3.3land Theorem 14.41 



Example 5.1. Let us consider the pure fractional equation dtu = g[u]. From 
e.g. [2H]; it follows that the solution of this equation is given by the convolution 
product u{x,t) = {K * uo)(a;,t), where K is the kernel of g. Using the properties 
of the kernel, it can be shown that this equation has a regularizing effect on the 
initial datum (see e.g. p'); this regularization appears clearly in our numerical 
experiments presented in Figure [l] 

Example 5.2. Let us consider the fractional transport equation dfU + d^u — g[u]. 
Our numerical results suggest that, as done by dtu + d^u = d^u, this equation regu- 
larizes and transports the initial datum. Our numerical experiments are presented 
in Figure [2] The numerical flux (3.3) has been used. 

Example 5.3. Let us consider the fractional Burgers' equation dtu + ud^u ~ g[u]. 
Our numerical experiments in Figure [s] confirm what has been shown by O : 
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-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



(a) T = 0.1 



(b) T = 0.2 



-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



(c) T = 0.1 



(d) T = 0.2 



Figure 2 . Initial data (piecewise linear) and solutions of the frac- 
tional transport equation (A = 0.5) with fc = and Ax = 1/160. 



this equation does not regularize the initial condition. Discontinuities in the initial 
datum can persist in the solution, and shocks can develop from smooth initial data. 
Figure [4] shows how the behavior of the solution changes with A: as A — >■ 0, our 
numerical solution approaches the solution of the pure Burgers' equation with a 
source, dtu + udxU = u; as A — )■ 1, our numerical solution approaches the smooth 
solution of the fractional Burgers' equation with A = 1 (see [IS]). Figure [s] clearly 
shows how a shock can develop and vanish in a finite time. Figure [6] shows how 
the accuracy improves with fc = 0, 1, 2. A third order Runge-Kutta (RK3) time 
discretization and slope limiters (cf. [M]) have been deployed in Figure |6] We 
have used the Lax-Friedrichs flux 

F(a, b) = i[/(a) + /(6) - c{b -a)], c = max{|/'(a)| : \a\ < ho||L-(K)}- 

Let us note that the above numerical flux does not fulfil assumption Al. However, 
this assumption can be replaced with a milder one: it is enough to ask F{a^ b) to 
be Lipschitz continuous on {{a,b) : \a\ < ||uo||l=°(r) and \b\ < ||uo||l=°(K)}- 

To give an idea about the speed of convergence of our experiments, we have 
computed their rate of convergence in Table [T] We have measured the error 

Eax.p ■= \\uAx{-,T) — {te(-, T)llip(R) 

{ue is the numerical solution which has been computed using Ax — 1/640), the 
relative error 

RAx,p := £'Aa;,p/||We(-, r)llLP(B), 
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-O.B -0,6 -0.4 -0,2 □ 0.2 0,4 0.6 O.B 



(a) uo{x) = -sgn(a:) 




■O.B -0.6 -0.4 -0.2 0.2 0.4 0.6 O.B 



(b) uq{x) = — arctan(15i')/90 




■O.B -0.6 -0.4 -0.2 0.2 0.4 0.6 O.B 



(c) uo{x) = sgn(x)l|3.|>i/4 + 4xl|^|<i/4 




-O.B -0.6 -0.4 -0.2 0.2 0.4 0.6 O.B 



(d) uo{x) = sin(27ra::) 



Figure 3. Initial data and solutions of the fractional Burgers' 
equation (A = 0.5) using fc = 0; T = 0.5 and Ax = 1/160. 



Table 1. fc = (left) as in Figure [s] (c) and k = 1 (right) as in 
Figure |6] (6). 



Ax 


Eax,1 


R-Ax,l 


"Ax,! 


Eax,2 


R'Ax,2 


aAx,2 


1/10 


0.1990 


0.1109 


0.5726 


0.4580 


0.3765 


1.0714 


1/20 


0.1338 


0.0746 


0.4711 


0.2180 


0.1792 


1.2024 


1/40 


0.0965 


0.0538 


0.3964 


0.0947 


0.0779 


1.1717 


1/80 


0.0734 


0.0409 


0.4399 


0.0421 


0.0346 


1.0881 


1/160 


0.0541 


0.0301 


0.7235 


0.0198 


0.0163 




1/320 


0.0327 


0.0183 











and the approximate rate of convergence 

aAx,p := (log-EAi^.p - l0g£'Aa;/2,p)/l0g2. 

We expected to see numerical convergence of order 1/2 for fc = and numerical 
convergence of order 3/2 for A; = 1 (i.e, high-order convergence). The values a^x,! 
roughly suggest 1/2 convergence while the values aAx,2 do not reach the expected 
rate 3/2. This could be due to our way or reducing the problem from a nonlocal to 
a local one (cf. Remark 5.1 1. 
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-O.B -0,6 -0.4 -0,2 □ 0.2 0,4 0.6 O.B 



(a) A = 0.1 




-O.B -0.6 -0.4 -0.2 0.2 0.4 0.6 O.B 



(c) A = 0.7 




-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



(b) A = 0.3 




-O.B -0.6 -0.4 -0.2 0.2 0.4 0.6 O.B 



(d) A = 0.99 



Figure 4. Initial data and solutions of the fractional Burgers' 
equation for different values of A using k ^ 0; T ~ 0.5, Ax = 1/200, 
and uo{x) — — arctan(15a;)/90. 



Appendix A. Technical lemmas 
Lemma A.l. Let (p e L^(M) n BV{M.). Then, there exists C > such that 

\\9M\LHm <cj f + < c,cMi-l;^ \^\^ 

JrJ\z\>o 1^1 

Proof. For all e > 0, 



\z\<eJm •'|z|<l 



\(p(x + z) — ip(x)\ , , „ \,, ,, f 1 , 

>' dxdz<2e-'yU.^^^ j-^dz. 

|z|>c JR Pi kl>l 1^1 

Set e = "*iiM to conclude. □ 
Lemma A.2. Let (f,(l) e L\M.) D BV{R). Then 

ipg[(t)\dx = / g[ip\(l)dx 

and, in particular, 

cx f f i^jz) - ipix))\ 
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-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



(a) T = 0.1 



(b) T = 0.7 



-O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 




■O.B -0.6 -0.4 -0.2 □ 0.2 0.4 0.6 O.B 



(c) T = 1.7 



(d) T = 2.9 



Figure 5 . Initial data (piecewise linear) and solutions of the frac- 
tional Burgers' equation (A ~ 0.5) at different times T using fc = 0; 
Ax = 1/200. 



Proof. By Lemma A.l and the fact that BV{R) C L°°(E), 



and, then, Fubini's theorem can be used to obtain 



< oo, 



ip{x)g[(f>{x)]dx 



1 



\z ~ x\ 



1 + A 



dzdx 



g[ip{x)](j){x)dx. 



□ 



Corollary A. 3. Lemma A. 2 holds true for all ip,(j) £ H^/'^{R). 



Proof. Lemma ] A. 2| holds true, in particular, for all 0„ step functions with com- 
pact support, 



(pnix)g[(f>n{x)]dx = / g[(pn{x)](j)n{x)dx 



(A.l) 



Let us choose, by density, (pm4>n — ^ (p, in i7^/^(M), and recall the following 
definition of the i/'^/^-norm (cf. [HI Chapter 6]): 



(A.2) 
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-O.B -0,6 -0.4 -0,2 □ 0.2 0,4 0.6 O.B -0,B -0,6 -0.4 -0,2 0.2 0,4 0.6 O.B 



(a) = (b) fc = 1 




-O.B -0,6 -0.4 -0.2 0.2 0.4 0.6 O.B -O.B -0.6 -0.4 -0.2 0.2 0.4 0.6 O.B 



(c) k = 2 (d) Solution computed using Ax = 1/640 

Figure 6. Initial data and solutions of the fractal Burgers' equa- 
tion at T = 1/10 using different values of fc = 0, 1, 2; Ax = 1/10, 
and uo{x) — sin(27ra;). 



Note that, using ( |A.2[ ) and ^ \ 

\\g[(pn] - 5MIIh-V2(r) = 



< 



{i+e)^^'['Pn{o ~ mM = ii^n - VWh 



X/2( 



since (1 + ^2-)-a/2^2A < (i + ^2^^/2 f-Qj. aU ^ e M (indeed, call = x, and muhiply 
both sides by (l+a;)^"'^/^ to get x^ < (l + a;)'^ which holds true for all x > 0). Thus, 
since g[ip n], g[ (l>n] siv'ljSM in i/"^/^(M) whenever ipn,4>n 'pA in 
equality ( |A.l ) holds true also in the limit n — > oo. □ 



-Lemma A.4. //(/)£ V^CM'^{M), then e H^I'^{M) and, for some constant c > 0, 



< 



Proof. Let us choose a function </> G C^ -^^(K), (t>{^) = X^igz Sp=o 
and let be the regular part of its derivative, 
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In this case we may define the quadratic variation of (p as 



First of all, we prove that both H^rlU^lR)) I'/'Iqv(r) are finite since (j) G L^{M.). 
Indeed, by orthogonality of the Legendre polynomials, 

Mhm^HY.^^^'i^^ andhence ^ ^p,^ ^ ^^^ll'^lli^(R)■ 
As a consequence 0^ € L^(R) since for each Legendre polynomial ipp^i, ^'Pp,i = 
apipp-i,i for some constant ap. Moreover, 

2 (fc + l)(2fc+l)„_^„2 

\v\Qvm) - 



IQV(R) ^ IIV^IIL^ 

while ¥>p,.(a;+) = (-I)p) 



since l^lgyfR) < 2 E*gz <^^(2;»^)+2 E»ez <^^K ) and (remember that = 1 



< (fc + 1) 2^ 2^ < II<^IIl2(m)- 

iez p=o 

Next, we prove that there exists a constant c > such that, for a.e. \z\ < 1, 
(A.3) [ [<j>{x + z) - c^{x)]^dx < c (\zMlvm + \^\^W\h(R)) ■ 

JR ^ ^ 

Note that 

r{i+l)\z\ 



I [(t){x + z) - (t){x)fdx ["^(^ + ^) - 

Jr iez-'^l^l 

= / ^[(^(.T + (z + l)|z|)-(/.(a; + i|z|)]2dx. 



By appropriately adding and subtracting the values 4'{xf), z e Z, the right-hand 
side of the above expression is less than or equal to 

(A.4) 3 r J2mxt) - H^7)?dx + 3 /"''' Y ['^(^P - ^i^j+i)^'dx, 
•'O iez •'O iGZ i=o 

where the Ji + 1 points x^ ~ Zq < . . . < Zj, — a;,^]^ lie inside the interval li (these 
points can vary from interval to interval depending on the value of Ax. E.g. if 
|z| ^ Ax, each interval li contains more than two points, while if \z\ 3> Ax, some 
intervals contain just the end-points xf = z'q and z}. = x~|_]^. We can control the 
first term in (A.4) thanks to the bound on the quadratic variation of (j) while, since 
inside each interval li the function </> is smooth, we can use the Taylor's formula to 
rewrite the second term as 

./i-i 

Y E [^^(4) - =^ 1^1 E E K(y;)]'(4+i - ^])' ^here z} < < z]^,. 

iez j=0 igZ j=0 

The right-hand side of the above inequality contains a Riemann sum approximation 
of the L^-norm of the function 4)'^ ^ C\ L?{E) and is therefore finite. Hence 
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inequality (A. 3 1 has been established, and we are now ready to conclude the proof. 
The seminorm 



kl>l 



[(j){x + z) — 4'{x)] 



1 + A 



-dxdz 



is finite since 



Ji < mihm 



\z\<l 



\z\>l 



[0(2; + z) - (/'(.t)]2 
|z|i+^ 



dz 



dxdz := Ji + J2 



and, thanks to (A. 3 1, 



J2 < c 



dz 

\z\<l \Z 



X \\K\\h(R) 



\z\<l 



dz 



< 00. 



□ 



Lemma A. 5. Let u G i7'^+^(M) anrf u 6e iis L?' -projection into , then there 
exists a constant c^, > such that 

II" ~ "llff^/2(R) — Cfe||w||ffA,+i(R)Aa;^'^^^ ^. 

Proof. Let us call v = u — and remember that (cf. [H Section 4.4]), for some 
constant c/c > and all intervals li = {iAx, {i + l)Aa;), 



||w||l2(/,) < Cfc||u||^fc+i(7.)Aa;''+\ 
||w||L~(/i) < Cfe||w||^fc+i(j_)A, 



fe+i 



a; 

||w||hi(/.) < Cfe||u||^fc+i(/^)Aa;''. 
First of all, let us bound from above the i/^/^-norm of v as 

|„||2 / |L,||2 , ^ f f [v{z)~v{x)]^ 



2 



/. /. „ ^^^^^ 



h Jh+l 

[v{z)-v{x)Y 

dzdx 

|2-x|>Aa: 



:= Ji + J2 + J3 + 



Note that 



Ji^Eii-"' 



We now prove the remaining Ji [i — 2,3,4) to be of order Ax^'"'+^ ^. First note 
that since v is smooth on each interval , the fundamental theorem of calculus 
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followed by Jensen's inequality yield 

l12 




li Jli 



\z-x[_ 
|1+A 

'+1 dzdx 

a" 



Thus h 



-JiJu k-^l' 

J Xi J Xi 

< MW.) / / ^ 

Jxi J-A: 



dsdzdx 



\z — x\ 
dsdx 

■Ax W 

ds 



= cAx^ '*'Z^i6z ll^llH-i(7i) < Cfc||u||^fc+i(^jAa;^'=+^ Next, we note that 



<2[ [ J^^dzdx + 2f f j^^.dzdx. 

estimate the first term on the right-hand side of the above 
ideas can be used for the second one. 

dzdx 



Let us show how to c 
inequality. Analogous ideas can 

A. A.+i \z-x\^+>^ 

< Ml 



\v 



J Xi J Xi^\ \^ 

rXi + 1 roo 

Hh+i) / / 

J Xi JXi+l—X 

TOO nx,. 



{z-xy+^ 

dsdx 



lz,~(/.+i) ^i+A 



S 

•Xi^i 



,1+A 



dx 



p dr r^'' dy 
cAa;i-^Eiezll^llioc(j,) < CkMl,+,^^^Ax^''+^-^. Finally, 

'|z-x|>Ax \Z-X\^+^ UJ\s\>Ax 

< 4||?;||i2(M) 



xy 



Thus /s : 



J'M.J\z-x\>Ax 



and I4 



= cAx-^Eiez IHl^ii) < c,\\u\\l,^,^^^Ax'''+^-\ 
Lemma A. 6. Let u gV'' (1 ap^i = fj g[u](fip^i and 

k 

fu{x) = ^^Op,i<^p,i(a;). 

constant c > 0. 



□ 



7« 



Then, ||7«||z,2(k) < c||u||z,2(r) for some 



24 



S. CIFANI, E.R. JAKOBSEN, AND K.H. KARLSEN 



Proof. Let us introduce the compactly supported function vm E O i^(M) 

k 

\i\<Mp=Q 

Note that, since ap^i = Jj g[u\ipp^i, 

\i\<Mp=0 \i\<Mp=0 •'-f' •'K 



By Lemma A. 4 the pairing J^g[u]vM is less than or equal to 



k 



l|w||//A/2(R)||t;A/|l//A/2(R) < c||u||l2(r)||wm||l2(r) < c||u||i,2(r) I ^ 

\|i|<A/p=0 

for some c > 0. Hence, J2\i\<M Sp=o '^p.i — c||M|li2(R) and, in the Hmit Af — oo 

fc 

iez p=o 



□ 



Appendix B. Proof of Proposition 14.1 
Since :— J^lj.{x)g[li.{x)]dx, Lemma 



A.2 



returns 



G]= f li^ix)g[li^{x)]dx = f li^{x)g[liAx)]dx^Gl 
Jr Jr 

Thus, by Lemma A.l X^jsz l^jl — /r Isfl-fi (^)] M^; < oo and, by symmetry. 

All diagonal elements are equal and negative. Indeed, 

li,ix + z)~li^{x) , _ f ({z 



where 

-\z\ z £ (-Ax, Aa;) 



—Ax otherwise. 

„1-A 



Thus, G:- = ~ca(/|^|^j i^'^-^^ + i|x|>i ji|T+xrf-z)Ax^ AU elements outside the 

- w 



diagonal are positive. Moreover, G^,X\ = G* for all e Z x Z since, if i 7^ j, 



. . l/.(x + z) 



Appendix C. Proof of Theorem 14.21 for the implicit-explicit method 
Let us consider the problem 
(C.l) V, - Atg{v), = hi, i e Z and /i e 1°°{Z) n /^(Z). 
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One can proceed as done by Droniou for nonlocal operators satisfying all the as- 
sumptions listed in [19] (cf. also [13] for a detailed proof for the operator g{-)) to 
prove the existence of a solution v G 1°°{Z) n of problem (C.l I. Moreover, 



(C.2) 
(C.3) 



inf hi < inf Vi < sup Vi < sup hi , 

ieZ iGZ jgx ieZ, 



Note that (C.2) ensures uniqueness for problem (C.l). Our plan is to rewrite 



the implicit-explicit method (4.1) in the form (C.l), and use (C.2)-(C.3| to prove 



Theorem 4.3 We start by rewriting (4.1) in "linearized" form, 



U. 



71 + 1 



where 



At ^ At Fiu^_„un-Fiu^,un 

Ax c/;;i - J7f ' Ax - up 



(the above coefficients are equal to zero when the denominators are equal to zero) 



By the CFL condition (4.31 and the Lipschitz regularity of F, it follows that 



ai,bi,l — ai — bi are bounded and positive. Thus, the implicit-explicit method (4.1 ) 
reduces to (C.l I if we choose Vi :— U^^"^ and hi 



a,C/«+i + (l-a,-fe,)C^r + fe^C^r-i- 

We are now ready to prove Theorem |4.2[ Items i and ii are easy consequences 



of (C.2) and (C.3). To prove item in, we call V^^ = U^^^ — C/f, and, by using the 
implicit-explicit method (4.1) in linearized form, we obtain 



n+l 



Ax^ ^ ^ Ax 

a^+lV:U + (1 - - h+i)V:' + b,Vr^,. 



Note that = E.ezGU^^^' = E.ezG^C/]!^ since Gjtl = Gj 

for all {i,j) £ Z X Z, and 



yn+l _ ^i T/n+1 

' Ax ^ 



G]V;+' = a,+,Vr+, + (1 - a, - b.,+,)Vr + b^V^.: 



which is of the form (C.l ) and, thus, Z^-contractive. This proves item in. The proof 



of item iv goes as the one for the fully explicit method (4.2) 
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Appendix D. Proof of Lemma [4l3] 



Note that A£^5[u, u] > by (4.8), and hence A^,s[u,u\ < A^j[u,u\ + A^.^fu 
Ii+ I2 + I3 + I4, where 



A / / Vu(y.s)iuix,t))ipt{x,y,t,s)dxdtdyds 
" iQt 

Vuiv.s) t))ipt{x, y, t, s)dxdtdyds, 
Quiy.s) t))(p^{x, y, t, sjdxdtdyds 

IQt JQt 

<lu{y,s){u{x, t))ipx{x, y, t, s)dxdtdyds, 



^s -"^ I I Vu(y,s)i'^i^^'t))9[u{x,t)]Lp{x,y,t,s)dxdtdyds 

'n'i(y,s){u{x, t))g[u{x, t)](p{x, y, t, s)dxdtdyds 



Qt JQi 
and 

h -^ / Vu{y,s)iu{x,0))^p{x,y,0,s)dxdyds 
Jqt Js. 

Vu{y,s) T))ifi{x, y, T, s)dxdyds 
Vu(y,s) {u{x, 0))(p{x, y, 0, s)dxdyds 

Vii(y,s) {u{x, T))ifi{x, y, T, s)dxdyds. 

As shown in [23, Theorem 3.11], Ii = I2 = while 
(D.l) h < c{e + S + Ax) - \\u{-,T) -u{-,T)\ 

We now prove that I3 < 0. Note that, since g[u] e L^{Qt), 



IQt 

f 

'Qt 



\v's.{y s){u{x,t))\\g[u{x,t)]\ip{x,y,t,s)dxdtdyds < 00, 

IQt JQt 

and we can change the order of integration to obtain 

/ v'u(y s)iHx,t))g[u{x,t)]ip{x,y,t,s)dxdtdyds 

JQt JQt 

+ / V'u(x,t)i^iy^s))9[uiy,s)](pix,y,t,s)dxdtdyds. 

JQt JQt 

Since r]'^{u) = -r]'^{u), 



/a = / / / sgn{u{x,t) -u{y,s))(p{x,y,t,s) 

JQt JQt J\z\>0 

{u{x + z,t)- u{y + z, s)) - {u{x, t) - u{y, s)) 
\u(x + z, i) — u{y + z, s)| — \u{x, t) — u{y, s 



dzdxdtdyds 



< 'Pix,y,t,s 

JQt JQt J \z\>0 



1^1 



dzdxdtdyds. 
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Let US rewrite the right-hand side of the above inequahty as a sum of two integrals, 
and use the change of variables (z, x, y) — > (— z, x + z,y + z) to obtain 



Lp{x + z,y + z,t, s) 

\z\>0 

\u{x,t)-u{y,s)\-\u{x + z,t)-u{y + z,s)\ 
. dzdxdtdyas 

^ I / f{x,y,t,s) 

Qt JQt J\z\>Q 

\u{x + z,t)- ujy + z, .^)| - Hx, t) - u{y, s) \ 



2 



1^1 

By adding up these terms we find that 

h<\ ( I I {fix + z,y + z,t,s) - f{x,y,t,s)) 

^ JQt JQt J\A>0 

\u{x,t)-u{y,s)\-\u{x + z,t)-u{y + z,s)\ 
dzdxdtdyds, 

and hence /s < since ip{x + z,y + z,t, s) = fix, y, t, s). To conclude, let us point 
out that the following result is needed in [231 Theorem 3.11] to prove (D.l|. 

Proposition D.l. Let u be a BV entropy solution of Jj. i^ . Then, there exists a 
constant c > such that t + S) — «(•, t)||^i(R) < c6. 

Proof. Let < a < 6 < T and f,] : M — > M be a smooth approximation of l[a^f,]. 
Let us call ip'{x,t) = (j){x)V^^ f^^{t), where e C^{R). Thus, 

uipl + f{u)(pl + ug[(p'^]dxdt — 

since u is a BV entropy solution of (1.1 1 and, so, a weak solution (cf. [2] for the 
definition of weak solution). The limit for e — > is, cf. [531 Theorem 7.10], 



and 



(t){x)[u{x, a) — u{x,b)]dx + / / f{u)(f>x + ug[(l)]dxdt = 



- M(-,a)||2,i(R) = sup / (t){x)[u{x,b) - u{x,a)]dx 

\<p\<l JR 

sup < - / / f{u)(j)x + ug[(j)]dxdt\ 



I0l<i 



ug[(l)]dxdt >. 



< c\uo\BV(M)ib- a) + sup < - / 

|0|<1 [ Ja 

To conclude the proof, the following estimate is needed: 

sup \ ^ / ug[(l)]dxdt > — sup \ — / (j)g[u]dxdt> < / [dxiif 

|0|<1 [ Ja JR J |0|<1 [ Ja Jr J Ja Jr 

< c{b-a), 

where Lemma I A . 2 1 and Lemma lA. II have been used. □ 



28 



S. CIFANI, E.R. JAKOBSEN, AND K.H. KARLSEN 



Appendix E. Acknowlegement 

We would like to thank the referees who did a very careful job reading this paper. 
Their indications and suggestions have not only contributed to make our exposition 
clearer, but have also helped us to correct and improve our original results. 

References 

[1] H. Abels, K. Kassmann. An analytic approach to purely nonlocal Bellman equations arising 
in models of stochastic control. J. Differential Equations, 236(l):29-56, 2007. 

[2] N. Alibaud. Entropy formulation for fractal conservation laws. J. Evol. Equ., 7(1):145-175, 
2007. 

[3] N. Alibaud, J. Droniou and J. Vovelle. Occurence and non-appearance of shocks in fractal 

Burgers equations. J. Hyperbolic Differ. Equ., 4(3):479-499, 2007. 
[4] A. Almendral, C. W. Oosterlee. Accurate evaluation of European and American options under 

the CGMY process. SIAM J. Set. Cornput., 29(1):9.'? 117, 2007. 
[5] P. Biler, T. Funaki, and W. A. Woyczyiiski. Fractal Burgers equations. J. Differential Equa- 
tions, 148(l):9-46, 1998. 
[6] P. Biler, G. Karch, and W. A. Woyczyriski. Asymptotics for conservation laws involving Levy 

diffusion generators. Studia Math., 148(2):171-192, 2001. 
[7] M. Bossy and B. Jourdain. Rate of convergence of a particle method for the solution of a ID 
viscous scalar conservation law in a bounded interval. Ann. Probab., 30(4): 1797-1832, 2002. 
[8] L. Brandolese and G. Karch. Far field asymptotics of solutions to convection equation with 

anomalous diffusion. J. Evol. Equ., 8(2):307-326, 2008. 
[9] S. C. Brenner and L. R. Scott. The mathematical theory of finite element methods Texts in 
Applied Mathematics, 15, Springer- Verlag (NY), 1994. 
[10] M. Briani, C. La Chioma, R. Natalini. Convergence of numerical schemes for viscosity solu- 
tions to integro-differential degenerate parabolic problems arising in financial theory. Numer. 
Math., 98(4):607-646, 2004. 
[11] M. Briani, R. Natalini. Asymptotic high-order schemes for integro-differential problems aris- 
ing in markets with jumps. Commun. Math. Sci., 4(l):81-96, 2006. 
[12] M. Briani, R. Natalini, G. Russo. Implicit-explicit numerical schemes for jump-diffusion pro- 
cesses. Calcolo, 44(l):33-57, 2007. 
[13] S. Cifani. PhD Thesis. In preparation. 

[14] B. Cockburn. An introduction to the Discontinuous Galerkin Method for convection- 
dominated problems. Lecture Notes in Math., 1697, Springer (Berlin), 1998. 

[15] R. Cont, V. Ekaterina. A finite difference scheme for option pricing in jump diffusion and 
exponential Lvy models. SIAM J. Numer. Anal, 43(4):1596-1626, 2005. 

[16] R. Cont, P. Tankov. Financial modelling with jump processes. Chapman & Hall/CRC Fi- 
nancial Mathematics Series, Chapman & Hall/CRC, Boca Raton (FL), 2004. 

[17] A. Dedner, C. Rohde. Numerical approximation of entropy solutions for hyperbolic integro- 
differential equations. Numer. Math., 97(3):441-471, 2004. 

[18] Y. d'Halluin, P. A. Forsyth, K. R. Vetzal. Robust numerical methods for contingent claims 
under jump diffusion processes. IMA J. Numer. Anal., 25(1):87— 112, 2005. 

[19] J. Droniou. A numerical method for fractal conservation laws. Math. Comp., 79:95-124, 2010. 

[20] J. Droniou, T. Gallouet and ,J. Vovelle. Global solution and smoothing effect for a non-local 
regularization of a hyperbolic equation. J. Evol. Equ., 4(3):479-499, 2003. 

[21] ,1. Droniou and C. Imbert. Fractal first order partial differential equations. Arch. Ration. 
Mech. Anal, 182(2):299-331, 2006. 

[22] G. B. Folland. Introduction to partial differential equations. Princeton University Press, 
Princeton (NJ), 1995. 

[23] H. Holden and N. H. Risebro. Front Tracking for Hyperbolic Conservation Laws. Applied 
Mathematical Sciences, 152, Springer, 2007. 

[24] G. Karch, C. Miao, and X. Xu. On convergence of solutions of fractal Burgers equation toward 
rarefaction waves. SIAM J. Math. Anal, 39(5):1536-1549, 2008. 

[25] A. Kiselev, F. Nazarov and R. Shterenberg. Blow up and regularity for fractal Burgers equa- 
tion. Dyn. Partial Differ., 5(3):211-240, 2008. 

[26] N. N. Kuznetsov. Accuracy of some approximate methods for computing the weak solutions 
of a first-order quasi-linear equation. USSR. Comput. Math. Phys., 16:105-119, 1976. 

[27] N. S. Landkof. Foundations of modern potential theory. Springer- Verlag (NY), 1972. 

[28] P. Levy. Calculul des Probabilites. 1925. 

[29] A. Matache, C. Schwab, T. P. Wihler. Fast numerical solution of parabolic integrodifferential 
equations with applications in finance. SIAM J. Sci. Comput., 27(2):369-393, 2005. 



DG METHOD FOR FRACTAL CONSERVATION LAWS 



29 



[30] M. Matalon. Intrinsic flame instabilities in premixed and nonpremixed combustion. Annu. 
Rev. Fluid Mech., 39:163-191, 2007. 

(Simonc Cifani), Department of Mathematic;s, Norwegian University of Science and 
Technology (NTNU), N-7491 Trondheim, Norway 
E-mail address: simone.cifaniamath.ntnu.no 
URL: http : //www . math . ntnu . no/~cif ani/ 

(Espen R. Jakobsen), Department of Mathematics, Norwegian University of Science 
AND Technology (NTNU), N-7491 Trondheim, Norway 
E-mail address: erj9math.ntnu.no 
URL: http://www.math.ntnu.no/~erj/ 

(Kenneth H. Karlsen), Centre of Mathematics for Applications (CMA), Department 
OF Mathematics, University of Oslo, P.O. Box 1053, Blindern, N-0316 Oslo, Norway 
E-mail address: kennethk8math.uio.no 
URL: http : //www . math . uio . no/"kennethk/ 



